Two dimensional electron gas near full polarization 
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We establish the consistency of the Fermi liquid description and find a relation between Fermi 
liquid constants for the two dimensional electron system near the point of full polarization due to a 
parallel magnetic field H. Our results enable us to predict connections between different thermody- 
namic properties of the system. In particular, we find that near the point of full polarization He, the 
thermodynamic compressibility of the system experiences a jump with the subleading {He — HY^^ 
dependence on the magnetic field. Also, the magnetization has a cusp with the dependence of the 
type {H - He) + {He - Hf^"^ at H < He- 

PACS numbers: 71. lO.Ay, 73.21. Fg 



I. INTRODUCTION 

The appearance of the new generation of high mobility 
heterostructm'es^ resulted in observations of interesting 
phenomena^^^ which brought new spotlight on the ef- 
fects of strong interactions in two-dimensional electron 
systems (2DES). For us this renews theoretical interest 
in the thermodynamical properties of clean Fermi liquids. 

Traditional theoretical description of 2DES is based 
upon the theories of weakly interacting Fermi gas^ and 
Landau Fermi liquid theory®'^. In the limit of very 
strong interactions (corresponding to low electron den- 
sities) such a description breaks down. However the sit- 
uation changes if the low density 2DES interacts with 
another 2DES with much higher density. In this paper 
we present an example where this coexistance naturally 
occurs as a result of a strong spin polarization of a 2DES 
due to a high in-plane magnetic field. We find a con- 
sistent Landau Fermi liquid description for this system 
(despite the fact that a naive estimate of the plasma pa- 
rameter for the minority spin component yields a for- 
mally large value). The remarkable feature of our result 
is that in the close vicinity of the spin polarized state the 
perturbative expansion is possible in terms of the inverse 
gas parameter of the low density subsystem. This enbles 
us to determine the functional form of the dependencies 
of the 2DES compressibility, magnetization, and specific 
heat on the small density of the minority electrons. 

To understand the relation of this problem to the Fermi 
liquid theory at zero field let us recall the basic structure 



of the quasiparticle interaction functional (we will not 
write the trivial long-range Coulomb interaction term) 



Hint = ^FPp^ 



(1) 
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where p is the charge density, S'^'J''^ denote components 
of the spin density, and F^, F'^ are the corresponding 
Fermi liquid parameters. The charge (singlet channel) 
and the spin (triplet) fluctuations are decoupled, and 
SU{2) symmetry of the system guarantees F^g = F'^Sa 



1/3- 



If the magnetic field is applied along, say, the x - di- 
rection (in the plane of the 2DES), the SU{2) symmetry 
is reduced to U{1) and one may write for the quadratic 
part of the energy 



"Hint — 



[S^f+F^PpS^ + ^ J2 ^"5". (2) 



This means that the system can no longer be described 
by two constants. Now the reduced symmetry allows for 
four independent parameters. 

Simplifications are possible, however, with the further 
increase of the magnetic field, because oscillations of the 
spin density components Sy,z become gapped (the gap 
equals to the Zeeman splitting Ez). Therefore, for the 
description of low lying excitations with the energy much 
smaller than Ez, the last term in Eq. (2) can be ignored. 
Introducing deviations of densities of majority (minority) 
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electrons ni(2) = p/2 ± 5^/2, one obtains the two-fluid 
model 



1 2 , ^22 2 , 

n, + — — n, + -Fi2nin2 



(3) 



characterized by three independent parameters. 

Equation (3) suggests two questions: (i) what is the 
lowest density of minority electrons for which it is appli- 
cable; and (ii) whether the three constants of the model 
are indeed independent. The ultimate goal of this paper 
is to show that (i) the Fermi liquid description is con- 
sistent for any density of the minority electrons and (ii) 
there is a relation between the Fermi liquid constants for 
the vanishing density of the minority electrons. The only 
requirement for this description to be valid is that the 
fully polarized electron system is a stable Fermi liquid. 
The remarkable feature of this result is that in the close 
vicinity of the spin polarized state the perturbative ex- 
pansion in terms of the inverse gas parameter TivfiI is 
possible (here vpi is the Fermi velocity of the minority 
electrons). 

The remainder of the paper is organized as follows. 
In Section II we give a phenomenological description 
of the system near full polarization, present the main 
results and predict connections between different ther- 
modynamic properties of the system. Section III con- 
tains the microscopic derivation (justification) of the an- 
nounced results, first in an intuitive, then in Subsec- 
tions IIIB and IIIC, in a more rigorous manner. 



II. PHENOMENOLOGY NEAR FULL 
POLARIZATION 



due to modification of the majority density. The last 
two terms characterize the shift of the energies due to 
the magnetic field, iJ, and = g^isH is the bare Zee- 
man splitting, with g being the bare (non-renormalized 
by electron-electron interaction) Lande gi-factor, and fiB 
being the Bohr magneton. The quantity E'^(N) corre- 
sponds to the value of the magnetic field above which 
the magnetization is independent of the field. In other 
words, this value limits from above the region of the field 
where the finite density of the minority electrons is still 
energetically profitable. 

In order to find the ground state of the whole system 
we have to minimize energy Eq. (4) with respect to the 
electron densities. Having in mind that the total electron 
density is fixed by an external gate, we note that the 
densities are coupled by the constraint 



6ni +n2 = SN, n2 > 0, 



(5) 



where 6N is the change in the total electron density with 
respect to the density threshold for the population of 
the minority subband controllable by the variation of the 
gate voltage. This yields either n2 = or 



+ 



dn2 
1 + Fn 



(6) 



(n2 -6N)-— (2n2 - 6N) . 



The crititical field E'^{N + 5N) is determined as the field 
at which n2 = solves Eq. (6). Then the first term 
vanishes (see below) and we obtain 



A. Structure of the theory 

Let us consider the system first at zero temperature. 
Because the total spin of the system commutes with the 
Hamiltonian, we can write the energy density of the sys- 
tem 5T-L in terms of the majority (rii) and minority elec- 
tron density (n2). Omitting the trivial term of the direct 
Coulomb interaction (we will work with fixed total den- 
sity 111 -\- n2 = N), one finds 

n = ^ t Snl + H2(n2, N) + —n26ni - 



(fei 



■ n2) - f TV. 



(4) 



The first term in Eq. (4) is the quadratic expansion of the 
ground state energy of the fully spin polarized electron 
system, n\ = N. It has the standard Fermi-liquid form 
with vi being the density of states (entering the slope 
of the specific heat). The second term H2 characterizes 
the energy of the minority electrons at fixed ni = N and 
the third term characterizes the change in this energy 



dEi 

dN 



1 + i^ii - F, 



12 



(7) 



Further progress requires knowledge of the function 
H2{n2,N). We find 



n-2 n 



H2{n2,N) = j dn' j dn" 



+ 



F22in",N) 



U2{n",N) vi{N) 



(8a) 

where the dependence of the density of states of the elec- 
tron density is given by 



+ 



V2{n,N) V2{N) 

+ Fii{N) - Fi2{N)Y 



+ - 



(8b) 



and the Fermi liquid constant for the minority electrons 
is 
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F22{n,N) = - [l + Fn{N) - 2F,2{N)] + 



(8c) 



+ 



4Vi^[l+Fn(jV)-fi2(^)]^ 



+ 



Equations (8) constitute the key point of this pa- 
per. They state, that properties of the minority elec- 
trons can be expressed in terms of the density of states 
i/2(-/V) = i^2{n2 — 0, A^) at the point of full polariza- 
tion (renormahzed by interaction with the majority elec- 
trons) and the Fermi liquid constants, Fn, F12 of the 
majority electrons. We will see below that the above 
equations impose certain connections between different 
observable quantities. It is interesting to notice that the 
relevant expansion parameter here is not the strength of 
the Coulomb interaction, e^, but rather its inverse power. 
This expansion is valid for n2a^ <C 1, where is the 
usual screening radius in two dimensions, (Ib — l/vie^- 

Postponing a rigorous derivation of Eqs. (8) until the 
next section, we discuss their physical meaning. Consider 
n2 — > and retain only the first line in Eq. (8c). Substi- 
tuting the result into Eqs. (8a) and (4) and keeping only 
quadratic terms we find 



ii , (i+i^ii) 

2V2 



2ui 



H = + ^' i ' (5nl - nl) + ^^n2 (n2 + Srii) . 

fx 



Let us now set = —71,2 (keeping electrical neutrality). 
The Hamiltonian takes the form 



n = 



"■2 

2V2 



which correponds to the compressibility of non - interact- 
ing minority electrons. This is not accidental; the major- 
ity electrons screen the Coulomb interaction at distances 
of the order of the screening radius. At small densities, 
however, the distance between minority electrons is much 
larger than this radius. Therefore, this screened interac- 
tion is seen by minority electrons only as contact interac- 
tion, i.e. the effect of this interaction vanishes because of 
the Pauli principle. This explains the origin of the first 
line in Eq. (8c). The second term in Eq. (8c) describes the 
effect of the finite interaction range. This effect clearly 
vanishes as the distance between minority electrons in- 
creases. Because the residual interaction is of the dipole 
type Usc{r) — e^a^/r^, its effect can be estimated as 
STi. ~ Uscir = 1/ s/ri2)n2 which immediately gives a ^/n2 
dependence to the Fermi liquid parameter F22 ■ As a mat- 
ter of fact, the same i/ni arises in all angular harmonics 
of F22 (and F12 as well; see the following section). This 
leads to renormalization of the effective mass and, there- 
fore, the density of states in Eq. (8b). The residual inter- 
action is weak and therefore the pcrturbative treatment 
of the minority electrons is legitimate. 



To complete the calculation of the ground state energy, 
we use Eqs. (8) in Eq. (6) and find with the help of Eq. (7) 



n2 = no 



Rs 



^Rsuy^ + ... 



(9) 



e'^v2 



V2 



dN 



no = iy2{N) [E'^{N + 5N) - E,] {El{N + 5N) - E,) 
Here 0{x) is the step function, and 



7 ^ 



32 
9^ 



0.639 . . . 



(10) 



is a numerical coefficient. It is noteworthy that the sub- 
leading term in this dependence is singular near E^ = E^. 
Substituting Eq. (9) into Eqs. (8) and the result into 
Eq. (4), we find for the ground state energy £ 



S = l±IlLSN--^-I± + ^-l^, (11) 
2vx 2 2v2 b V2 ^ ' 

where we omitted terms linear in 8N (apart from the 
term proportional to E^) as they will not contribute to 
the observable quantities discussed in the following sub- 
section. 



B. Experimental consequences. 

In this subsection we apply the above ideas to establish 
relations between different thermodynamic properties of 
the system near the point of full polarization. 

The leading contribution to the specific heat of the 
two-liquid system is simply the sum of the quasiparticle 
specific heats of the two species of electrons. 



dCy 
dT 



+ 



+ 



(12) 



S-f^Rliy2SE, 



^{SE,): 



SE, = EliN) - E, 



This gives the operational definition of the density of 
states, even though the actual measurement of the spe- 
cific heat in two dimensions is technically difficult. 

Next, Eq. (7) allows one to find a certain combination 
of the Fermi-liquid constants 1 + Fn — F12 from the mea- 
surement of the critical magnetic field as a function of 
electron density. 

Further information about the Fermi liquid constants 
can be obtained from studying the thermodynamic com- 
pressibility K = {d'^£/d'^N)~^, where £ is the ground 
state energy of the system. Differentiating Eq. (11) yields 
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1 ^ 1 + Fn 

K Vi 



37 



^jRiu25EAe{5E,). (13) 



Wc sec that the measurement of the jump in the com- 
pressibiUty gives the value of the parameter Rg which to- 
gether with Eq. (9) and measurements of dNE^{N) yields 
the value of V2{N) without specific heat measurements. 
After that the coefficient in front of subleading square- 
root singularity docs not contain any fitting parameters. 

Finally, we calculate the magnetizetion M — —dnS as 
the function of magnetic field H. Differentiating Eq. (11) 
at SN = with the help of Eq. (9), we obtain 



M = gijLB 



N 

— - V25E^ 



{l-^.jR^,U25E,)e{5E,) 



(14) 

It is important to emphasize that after the compress- 
ibility and M at Ez > E'^ are measured, the formula for 
magnetization at Ez < E^ will not have any adjustable 
parameters. It is also worth noticing that the sub- leading 
dependence has a square root singularity, similar to that 
in Eq. (13). 

The predicted dependencies of magnetization and com- 
pressibility are plotted on Fig. 1. 



1+Fl 
Vl 



1+F 
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Vl 



1 

K 



A 



A 



dM 



g|.lBV2 



presented after Eq. (8c) and it will also be confirmed by 
the microscopic calculation. 

The route we are taking in this section is the following. 
First, we express our physical arguments in the language 
of the perturbation theory, i.e. identify the set of dia- 
grams leading to Eqs. (8). The rest of the derivation aims 
at showing that (i) these are the only diagrams that pro- 
duce a combination of the constant (n2 = term) and the 
term with y/n2 singularity, (ii) all other diagrams result 
in contributions of the order of {n2/N)lnn2/N; (iii) the 
Fermi liquid description of the minirity electrons is jus- 
tified. This material will be structured into subsections 
III A, IIIB and IIIC. The rather cumbersome content of 
these subsections will not directly contribute to the final 
results and might be skipped by a pragmatic reader. 



- 



A. Perturbation theory 



The goal of the microscopic consideration presented in 
this section is to prove the main assumption of the phe- 
nomenological treatment, namely Eqs. (8). We start by 
showing how one can arrive at Eqs. (8) using a simple- 
minded perturbative approach. The reason to do this 
is to build up physical intuition, identify the group of 
diagrams which gives the dominant contribution to the 
minority electron interaction, and to clarify the assump- 
tions, which one needs to make in order to justify this 
treatment. In the following two subsections we shall 
prove the validity of these assumptions and provide a 
more rigorous treatment of the problem. 

The main physical idea of the following picture is that 
close to the point of full polarization only the majority 
electrons participate in screening. Therefore, as a build- 
ing block for the perturbation theory we are going to 
use the dynamically screened (by the majority electrons) 
Coulomb interaction, which propagator has a diagram- 
matic representation of Fig. 2. The corresponding ana- 
lytic expression is 



FIG. 1: Predicted dependences of compressibility (dashed- 
dotted line) and susceptibility (solid line). 



III. MICROSCOPIC DERIVATION 

The purpose of this section is to develop the micro- 
scopic description leading to Eqs. (8), which is needed 
to justify the Fermi liquid description of the minority 
electrons and to calculate the coefficient in the second 
term of Eqs. (8b) and (8c). The form of the first term 
in Eq. (8c) follows already from the physical argument 



where Vo{q) = 27re^/|g| is the bare Coulomb interaction 
and Hi {LJ,q) is the polarization operator of the major- 
ity system. It is defined as the part of the density - 
density correlation function irreducible with respect to 
one Coulomb line. For small momentum and frequency 
transfers, q,u/vFi <C Pfi, it has the usual Fermi-liquid 
form*" 
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Hi (w, q) = J "^^^1 (w, q; ni, n2) ; (16) 
III (i^, q; ni, n2) = Uf (w, q; rii) 27r5(0i - ^2) - 



d4>3 
2n 



nf (a;,q; rii; 



nf (a;,q;n) 



vivpin ■ q 



vpi^ • q — w — iOsgnw 



Hi (w,q;n3,n2) ; 



(17) 



Here Fn (11^112) is the majority Fermi liquid parameter 
(taken at zero minority density n2 = 0), {uj,q\n) 
characterizes the hnear response of non-interacting quasi- 
particles, and v\ is the density of states of majority quasi- 
particles. The latter quantity enters into the specific heat 
of the spin polarized (n2 — 0) system. Equation (16) 
takes into account all possible contributions singular as 
a function of vpiq/uj and neglects all the contributions 
which are regular functions of the parameters {q/pFiY 
and Lojtpi. Unit vectors = (cos 0^, sin c/ii) characterize 
the direction of motion of the quasiparticle. 
Now we discuss the origin of Er^s. (8). 



<Cni)> 



FIG. 2: Dyson equation for the dynamical screening by ma- 
jority electrons. 



1. Interaction of minority electrons 

To justify the relation Eq. (8c) we need to dc;scribe 
the interaction of minority electrons in terms of the pa- 
rameters of majority electrons. Since interaction between 
minority electrons is characterized by the energy transfer 
w ~ qvF-2 "C qvF\, we can use the static approximation 
for F(w, q). Moreover, at wavectors smaller than the in- 
verse screening radius of majority electrons i^iVo((7) ^ 1. 
In this case Eq. (15) becomes (see subsection IIIB for 
further discussions) 



V{u:,q) 



1 + F^°^ - 



■ (1^ 



Here we use angular harmonics of the Fermi liquid func- 
tions 



(0) = E' 



(m) 



(19) 



for i,j = 1,2. The zeroth angular harmonics, F^j \ that 
appear in Eq. (18) correspond to the constants used in 
the previous sections. 



FIG. 3: Leading contribution to F22. 



The lowest order contribution to the minority electron 
interaction (that determines the Fermi liquid constant 
F22) is given by the diagram in Fig. 3. Using Eq. (18) to 
evaluate this contribution we find 



Fig. 3. = - 



1 + Fff 
1^1 



+ 



z^?^o(|pi -P2I 



(20) 



This contribution is proportional to |pi — P2I ~ Pf2 oc 
^Jn2 which produces the singular density dependence in 
Eq. (8c). 

We now show that higher order diagrams built out of 
the same ingredients as the simplest diagram in Fig. 3 
depend on at least the second power of |pi — P2I and 
thus result only in regular contributions to ^22- 

Indeed, summation of the ten second order diagrams 
in Fig. 4 yields zero whenever one of the dashed lines is 
substituted with a constant. Therefore, the constant part 
of the potential Eq. (18) can be omitted entirely from the 
second order perturbation theory and one obtains 



Fig. 4. (a) 



(g) oc 



1 + F 



(0) 



Z^fK)(Pl2) 



OC (Pl2)^; 



P12 = iPl - P2|- 



This cancellation is not accidental and in fact is due to 
the Fermi statistics of the minority electrons. All higher 
order terms are canceled in the same manner. 

Unfortunately, this is still not the whole story. Major- 
ity electrons affect minority electrons not only through 
the density-density interaction but also through the 
renormalization of the spectrum of minority electrons 
(the simplest analogy here is the polaronic shift of the 
bottom of the minority band). This renormalization de- 
pends on the distribution function of majority electrons 
and, therefore, generates the Fermi liquid function F12. 
The lowest order diagrams for this parameter are shown 
in Fig. 5. Precisely the same diagrams enter into the 
two particle irreducible vertex (that contributes to -F22) 
in Fig. 6. Therefore, this so far neglected contribution 
to the minority interaction can be expressed entirely in 
terms of F\i (by means of direct comparison with dia- 
grams in Fig. 5). We find 
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(a) 







(d) 



(b) 



(e) 

1 <^ 







-> ' ^ 



(g) 



>-o-< >■ 

(h) (i) 



-^J — 1 — ^ 



> 1 ^ 



(c) 

-< I <r- 



(f) 



(j) 



where i = 1,2. Summing up contributions of Eq. (21) 
and (20) we obtain the angle-dependent Fermi hquid pa- 
rameter F22 



F22{0) 



1+F 



(0) 



11 



[l + ^ff-fl2 

viVo{2pF2 sin ' 



(0)i2 



(22) 

The zeroth angular harmonic of Eq. (22) gives precisely 
Eq. (8c), if we recall the Landau theorem that relation 
1^2 = p'f2/^''^ is '^ot changed by interaction in any order 
of the perturbation theory. 



FIG. 5: Lowest order diagrams for F12. 



FIG. 4: Connected second-order minority diagrams (all solid 
lines denote Green's functions of minority electrons). 



2. Density of states of minority electrons 



Fig. 6. (a) 
Fig. 6.(6) : 



%^n,(o,pi.) = 



-la+i^i^) 



(21a) 



rp(0)i2 

^^ni(o,pi2)V(o,pi2) (21b) 



(0)12 



+ 



z/i(l + Fi(°)) z^iV^,(pi2) 
Fig. 6.(c) = Fig. 6.{d) = Fi^^U,{0,p,2)V{0,pi2) 



F. 



(0) 



J^i 



1 + fff 

J^lVb(pi2) 



(21c) 



The reason that only the zeroth harmonics of the Fermi 
liquid parameters appear in Eq. (21) is that we assume: 

(i) fI°'^ « const + 0{q^/pl-^) (where q is the transmitted 
wavevector); (ii) F^l^ ~ 0{q/pFi)- The dependence of q 
on the scale of the order of ppi can be neglected because 
it would generate smallness of the order of n2/N. 

Finally, to obtain the Fermi Uquid parameter F22 
we combine the two-particle irreducible minority vertex 
functions discussed above: 



Z^l 



= Fig. 3 + Fig. 6(a) + • • • + (d) 



\Pj[j=pf2; 
P1P2 = 0, 



In order to determine the DoS (or the effective mass) of 
minority quasiparticles. let us recall Galilean invariance, 
which results in the following two-liquid variant of the 
usual Ward identity 



1 _ 1 1 r de (pFi 



m2(n) 



m 



mi J 2tt {pf2 



Fi2i9)+F22ie)^ 00^(9), 
. (23) 

where m is the bare (band) electron mass, and mi is the 
majority quasiparticle mass. At zero minority density 
F\2 renormalizes the mass by the amount of order one 
[since the large factor pfi/pf2 cancels exactly due to the 
assumptions that we made deriving Eqs. (21); see the 
assumption (ii) in the paragraph followingEqs. (21)]: 



1 



m2(0) 



1 

m 



1 



PFl 



lim 

mi PF2^0 PF2 



de 
2^ 



Fi2ie)cos{e). (24) 



After the cancellation of the prefactor the remaining de- 
pendence of ^12' on 712 is analytic. Hence, this term does 
not produce any singular dependence of the minority DoS 
on 712- The square root dependence is caused entirely by 
the interaction between minority electrons, i.e. F22. By 
plugging in Eq. (22) to Eq. (23) and using Eq. (24) one 
immediately retrieves the DoS Eq. (8b). 



3. Underlying assumptions 

The line of argument presented so far relies on several 
assumptions which require further justification. These 
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(a) 




defined subclass of minority 2-particle vertices, (provid- 
ing us with the nonanalytic dependence for F22), all n- 
point functions are smooth (i.e. Taylor-expandable) as 
a function of external momenta. Second, we will treat 
n2 <C A'' as a perturbation and show that this may intro- 
duce only corrections linear in the small parameter n2/N. 
We will use the vertex functions and the gauge invariance 
of the theory to calculate the Fermi liquid constant F22 
in terms of the Fermi-liquid parameters of majority spin 
and obtain Eq. (22). Finally, we will justify the calcula- 
tion of the minority mass in more detail. 



> — I I > 

(b) 

< 1 1 < -< 1 <r- 




> 1 > > ' ' >— 

(c) (d) 

FIG. 6: Contributions to F22 from the renormalization of the 
minority spectrum. 



include: (i) the momentum dependence of F12 was as- 
sumed to have the specific form (see the text following 
Eqs. (21); (ii) the screened interaction V{LU,q) was only 
considered in the limit of small frequencies, based on the 
intuitive assumption lu ~ vf2<1 ^ vpiq [see text preced- 
ing Eq. (18)]; (iii) the consideration of minority interac- 
tion was limited to certain class of diagrams, see above. 

The assumption (i) allowed for the explicit result 
Eq. (22) that followed from the evaluation of the dia- 
grams in Figs. 3 and 6. The essence of the assumption 
(iii) is that no other diagram contributes to the singular 
dependence of F22 ■ Partially this was illustrated by con- 
sidering diagrams in Fig. 4, however, one could imagine 
more complicated diagrams involving majority electrons. 
Moreover, the cancellation of diagrams in Fig. 4 relied 
heavily on the assumption (ii). Thus, in order to rigor- 
ously prove our conjecture Eqs. (8) we need to justify 
the above assumptions. Although proper consideration 
of these issues will not change the final results, we in- 
clude the following subsections in order to complete the 
derivation. 

Our strategy will be the following. First, we will set 
the number of minority electrons 712 to zero inside dia- 
grams in consideration, and discuss the analytic and scal- 
ing properties of the self-energy and n-point irreducible 
vertex functions. We will see that apart from a well- 



Be Completely polarized system 

In this subsection we discuss the properties of the irre- 
ducible vertex functions which we shall use in the follow- 
ing subsection to calculate the Fermi liquid parameters. 

1. Green functions 

We start by defining zero-temperature, real time Green 
functions of minority {j = 2) or majority {j = 1) elec- 
trons: 

(T,^,-(Ri)V'i(R2)) (25) 

where all the operators are taken in the Heisenberg rep- 
resentation and averaging is performed over the ground 
state of the system. To shorten the notation we use here- 
inafter the (1 + 2)-dimensional vectors R = {t,r) and 
P = (e, p) with the scalar product PR = et — pr. There 
are no Green functions mixing the electron species be- 
cause the electron spin component along the magnetic 
field is conserved. 

Let us set the number of minority electrons n2 to zero. 
For majority electrons we will need only the linearized 
spectrum near the energy shell: 

Gi{e,p) = (26) 

\vfiPfi vfiPfiJ 

where ^1 = Ufi (|p| — pfi) is the distance from the Fermi 
surface, Zi is the quasiparticle weight, vfi is the Fermi 
velocity renormalized by interaction, andp^i is the Fermi 
momentum. The relation [pfi]"^ = 47rni (ni being the 
density of majority electrons) is not affected by interac- 
tion due to the conservation of the number of states and 
the spin conservation (Landau theorem). The remainder 
of the self-energy for the majority electrons possesses the 
following property: 
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Si(a;,t/) = ReEi(a;,2/) +zsgna; |ImSi(a;,y)| . (27) 

The leading dependence of the self-energy in two dimes- 

nions is T,i{x,x) ~ «|a;|.xln(l/|.x|). 

As usual in Fermi liquid theory, the leading non- 
analytic dependences of vertex functions orginate from 
the overlap of poles of two Green functions with close 
momenta. In this case we will use the standard represen- 
tation 



(28) 



nr(c^,g;n) 



1^1 



+ ^(P,Q), 



where 2+1 vector Q = (w, q) is small in comparison 
with the Fermi momentum, n = p/|p| is the unit vector 
along the momentum p, and </?(P, Q) is a smooth func- 
tion with well defined limit at Q — > 0. The quasiparticle 
polarization operator was defined in Eq. (17). 

The minority electrons are described in a similar man- 
ner with the exception that now the spectrum cannot be 
linearized: 



G2(e,p) 



e-f--S,(e,p) 



(29) 



where m is the bare mass of the electron and S2 is the 
self-energy of minority electrons. Because n2 = 0, condi- 
tion / (iee*'^'G2(e, p) = must be satisfied and therefore 
G2{e) is an analytic function of e at Ime > 0. The self- 
energy has the effect of renormalizing the residue, the 
mass and the chemical potential: 



-S2(e,p) = iO + Z^^fX2 + (^2"^ - l)e 



(30) 



1 

m 



2 



■^2 



Here the parameter Z2 is the quasiparticle weight for the 
minority electrons. It has the physical meaning of the 
overlap of the initial wave-function of the majority elec- 
trons with the wave- function of these electrons after they 
screen the potential of an introduced minority electron. 
The chemical potential ^2 is shifted with respect to its 
bare zero value by the interaction with the majority elec- 
trons. This effect is analogous to the polaronic shift of 
the bottom of the band. The same polaronic effect in- 
troduces the renormalization of the electron mass, m^- 
There are two important points worth mentioning here: 
(i) E2 cannot introduce linear in p corrections to the 
spectrum, and (ii) the sign of TO2 is unknown, we will 
assume that it is renormalized to a positive value. 



The parameters Z2, H2, and m2 (recall that we are 
discussing the case of 712 = 0) arc determined by the 
integration over large momenta of the order of pFi. That 
is why they can not be calculated from the first principles 
and we treat them as input parameters of the theory. If 
the interaction in the majority sector is weak, then the 
calculation of /X2, and rn2 is possible. The imaginary 
part of the remainder of the retarded self-energy can be 
presented in the form (see Appendix A) 



-ImS2 (e;p) 



'2e 



Pfi 



2m2\e\ 



\PfiJ 



(31) 

where f{x) is a dimensionless function with properties 



2/3 X 



All 



(32) 



The above self-energy describes in particular the finite 
lifetime of the minority electrons with respect to the 
emission of the electron-hole pairs in the majority liq- 
uid. The rate of this decay is proportional to e^/^. This 
is different from the usual In e for the two dimensional 
Fermi liquid because there is no Fermi surface for the 
minority electrons formed yet. However, the quasiparti- 
cles are still well defined even in this case provided that 
e <S VfWfx- The form of Eq. (31) follows from simple 
dimensional analysis of corresponding diagrams which is 
elaborated upon later in this section, however the di- 
mensionless function J{x) can be obtained only by direct 
calculation, sec Appendix A. 

As we already mentioned, the Green's function for the 
minority electrons at n2 = is an analytic function of e 
at the upper semiplane Ime > 0. Therefore, the contribu- 
tion of two close poles is not dangerous and the singular 
part of the type of Eq. (28) does not arise. 

Concluding this subsection, we emphasize that we have 
assumed that the curvature of the spectrum at A; = 
for the minority electrons is positive. All the further 
scheme is based on this assumption, which we cannot 
justify for arbitrary interaction strength. We will not 
speculate on the alternative scenario in this paper. For 
more information on the minority Green's function and 
self-energy we refer the reader to Appendix A. 



2. Vertex functions - general definitions 

To characterize interaction of the minority electrons 
with each other as well as with the majority electrons we 
will need 2n-point vertex functions, which we denote by 
■ (PJ«*, . . . P7*; P^", . . . Pj^"), and define as 
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-i{2TTf5 (pr* - pr) j Aii'.'.i: (prs • ■ • pr*; pr. • • ■ pd = j ^^R-r* • • • rfRr*'^Rr ■ • • ^^Rjr (33) 
X exp l^z (pr'Rr' - prRi")) {T^^n (Rr*) ■ • • V',. (Rr*)V'l (Rr ) • • • v^l (rd)"'"^- 

The external legs on each diagram are amputated. 

Because of spin conservation, the number of incoming legs with j = 1(2) equals to the number of outgoing legs 
with j = 1(2). Because of the Fermi statistics the vertex function is antisymmetric with respect to the permutations 
of outgoing legs 



_^il...in 
11. ..,11, 



ll,...,l„ 



Tyout 

' ik ' 



pr 



/'■pout "DOU 

[^1 ^"'^^il 



Tyout 



■pout -pout. Tiin Tiin\ . 



(34) 



and have the same property for the ingoing ones. 

Our aim is to identify the relation of the Fermi liquid 
constants with the vertex functions in the theory. To do 
that we follow the standard procedure of the Fermi liquid 
theory and explicitly separate those contributions to ver- 
tex functions that contain possible singularities. There 
arc two sources of singularities: (i) overlaps of poles of 
two Green's functions, see Eq. (28), and (ii) the Coulomb 
propagator. Exact vertex functions Eq. (33) can be built 
using n'"^ and Vd as the basic building blocks in addition 
to the nonsingular part of A: 



irreducible 



(35) 



Here "irreducible" means that V comprises all the dia- 
grams that can not be cut by one Coulomb line or two 
majority Green's functions. More precisely, in each di- 
agram reducible in two majority electrons we substitute 
only the smooth part of the product of two corresponding 
Green's functions 



where P = (e, p). Notice that due to the irreducibility 
definition, see Eq. (36) , the question of the order of limits 
is resolved automatically. 



A = 



■X 



X 



+ J (^^m^ 




Gi P 



Gi P 



9l 

2 



<^(P,Q), 



(36) 



see Eq. (28). Irreducible vertices (35) obey the antisym- 
metry relation (34). 

The total vertex function can be quite easily found 
from the irreducible ones. The corresponding relation 
for the 4-point vertex function is presented on Fig. 7. In 
this scheme the Fermi liquid parameters are going to be 
determined by the vertex W . 

The 3-point irreducible (in the same sense as F) vertex 
function Bi(P, Q) satisfies the Ward identity 



B,(P,0) = 1 



de ' 



(37) 



P-Q/2 



P + Q/2 



8(e)8®nf 



FIG. 7: Relation of irreducible (F, B) and reducible (A) ver- 
tices. 

At n2 = any closed loop for the minority electrons 
vanishes, so majority electrons obey the standard Fermi 
liquid description, which does not depend on the values 
of and rig. Calculation of vertices involving minority 
electrons, however, requires knowledge of the irreducible 
vertex functions and V"^; we will need their values at 
external momenta much smaller than ppi- We intend to 
prove the smallness of and determine the dependence 
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of on external momenta in this region. The proof is 
based on the dimensional analysis of each order of the 
perturbation theory for the minority vertices. 

In building further perturbation theory for the finite 
density of minority electrons, we will use the screened 
interaction (15) as the basic interaction propagator, be- 
cause it already contains all the singularities (28) of the 
theory. In particular, V{u! = 0, q) is always finite and 
short range, unlike the bare interaction. We will see later 
that the relevant contributions come from u> <C vpiQ', in 
this region, we can easily solve Eq. (16), and obtain from 
Eq. (15) 



V{q,u;) « - 
1^1 



1 + f(?) - 



i\uj\ 
qvFi 



+ 



_ 1 
The wavevector 



qvF\ 



+ ... 



(38) 



Fig. 8). We notice that there are two scales in the prob- 
lem. The first, "ultraviolet" scale is determined by the 
majority electrons, i.e. the wavevectors of the order of 
Pfx and the energies of the order of vfxPfx. The sec- 
ond, "infrared" scale is determined by the momenta and 
energies of the external legs, i.e. the scale of the integra- 
tion over the momentum and energy is qiven by Q and 
Q'^ji^m-i) respectively. Statement (40) is the obvious 
consequence of the perturbation theory for infrared di- 
agrams. Indeed, the lowest order nonvanishing diagram 
for the 2n-point function contains n interaction lines and 
n electron minority Green hmctions, see Fig. 8. Calculat- 
ing a diagram in this regime one can use approximation 
(38) for the interaction potential. The constant part of 
the potential (corresponding to contact interation) can- 
cels immediately from the whole theory (i.e. from any 
vertex function A) when being substituted in any inter- 
action line, see e.g. Fig. 4., due to the antisymmetry 
relation (34), which is a simple manifestation of the fact 
that spinless fermions are not sensitive to the contact 
interaction. Only does the remaining part of the interac- 
tion 



: VFl 



e 



(39) 



characterizes the screening of the Coulomb potential by 
majority electrons. For reasonable interaction strength 
qs is not that different from "ppi- That is why we will not 
write ratio qs/pFi in the subsequent estimates unless it 
is necessary for a quantitative analysis. 



3. Vertex functions V for minority electrons 

Let us start with the vertex function involving ingo- 
ing and outgoing legs for minority electrons. We intend 
to show that the form of the potential (38) and the an- 
tisymmetry relation (34), guarantees that for the small 
external energies and momenta, pi ^ Pfi, ^ vfiPfi, 
the 2n point function has the following structure: 



^2,...,2/ 



■ 2,...,2\^l 1 • • • I '^n ) 

Z^2(PFl)"0"-4 

n 

E([pr]' + [prf) 

.i=l 



Q = 



1- 

1/2 



.PFlJ 



(40) 



where 7*^"^ is a finite dimensionless function, obeying the 
antisymmetricity relation following from Eq. (34). 

Relation (40) can be shown as following. Consider any 
order of the perturbation theory (lowest non-vanishing 
diagrams for 4— and 6— point vertices are shown in 



6V{q,iv) = 



qs qvFi 



contribute to the final answer. Then any infrared integral 
like 



J„ = y d'pde[5V{p,e)r[G2{p,i 



[we omit external momenta here] can be made dimen- 
sionless by expressing all momenta in units of Q, and all 
energies in units of Q^/2m2- This way, we obtain 



2 Q' 



-V 



/ dimensionless\ 
2m^ \ 2m^ J \m2PFi J ^\ function J ' 



which clearly has the form of Eq. (40). Inclusion of addi- 
tional m interaction lines bearing small momenta into the; 
tree level diagrams in Fig. 8, provides additional small- 
ness (Q/pFi)™. 

This procedure of finding the scaling form of the mi- 
nority vertex function F relies on assumption that the 
integrals are determined by the small momentum region. 
To justify this assumption, let us show that the contri- 
bution from the ultraviolet is always small. Indeed, let 
us separate the contribution into F from Eq. (40) where 
all the integrals are determined by the ultraviolet parts. 
Then we can introduce the momentum scale k* = Appi, 
where A is a numerical coefficient smaller than 1, and re- 
strict integration over momenta by > fc* and over the 
energy by |e| > VFik*, and call this contribution F. Be- 
cause the integrals are restricted to the high momentum 
region, F is an analytic function of its external momenta 
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FIG. 8: Lowest order contributions to 4- and 6-point minority 
vertex functions F. Dashed lines denote the non-constant part 
of the screened Coulomb interaction SV, internal lines are 
minority propagators G2. Black triangles denote the vertex 
B, see Fig. 7. 



and energies and can be expanded in Taylor series. Be- 
cause of the antisymmetricity constraint (34), the first 
nonvanishing term has the form (put all external ener- 
gies to for the sake of simplicity): 



p2,. . .,2 /-ryOUt 



r 



(P 



pm 

1 ) ■ ■ ■ I -"^ n 



r i2n-4 

[Pfi\ Vi 



(41) 



[k 



•\n{n—'\.) 



Because, by construction, k* is different from ppi only 

by a numerical factor, this estimate is smaller than the 

2 ^ 

result of Eq. (40) by a factor of {Q/pfiT ~ , for n > 2. 

The case of the 4-point vertex is special - ultraviolet 
and infrared estimates have the same powers of p^i in 
the denominator. It means that this vertex function is 
uniformly contributed by all energy scales. Therefore, 
function 7^^^ , may contain logarithmic dependence of the 
high energy scale ppi- Indeed, the direct calculation of 
this particular vertex shown in Appendix B gives 



^2,2 
■ 2,2 



out T^out.-nin ■Tyin\ 

^2 1^1 '^2 j — 



1 

VI [p. 



(42) 



71 m 1- 72 + C 

Vf\ \Pfi 



where the dimensionless functions 7, describe the de- 
pendence on the external momenta and energies: 7^ = 

({^} ' {~^^})' From the dimensional analysis 
above, it follows that all the leading graphs should have 
one infrared loop similar to the tree level diagrams. All 
other loops must be ultraviolet their role is the "dress- 
ing" of the vertices of the tree-level diagrams [black tri- 
angles in Fig. 8]. This dressing changes the numerical 
coefficient in the final expressions but does not affect the 
analytic structure of Eq. (40). 



4- Vertex functions F involving minority and majority 

electrons 



The next object to consider is the vertex function in- 
volving both minority and majority spins. We start from 
the simplest vertex with two minority and two major- 
ity legs. This object characterizes the correction to the 
simple Coulomb interaction between minority and ma- 
jority electrons, the lowest order diagrams contributing 
to this vertex are shown on Fig. 9. We will be interested 
in the behavior of this vertex when the energy and mo- 
mentum transfers are small in comparison with the Fermi 
energy and Fermi momentum of majority electrons, but 
are arbitrary in comparison with energy and momentum 
of minority electrons. Furthermore, we are interested in 
the situation, where the majority legs are nearly on-shell. 
We intend to show, that this vertex has the form 



u,Tll{Pr, Pr*; Pi", Pl") = lif + ^^V^i^ (43) 

[PFl] 



Pl':2= (ei,2± 2' Pl,2± 2 



P2 

PF\ 



T>OUt 

1,2 



ei,2 T ^, Pl,2 T ^ ) , i\=VF\{\V\\-VF\) 



where 7^2^ and i]"^^ are finite numerical coefficients. 

To understand the relation (43), we notice that it is 
equivalent to the statement that can be expanded 
in a Taylor series as a function of momenta and energy 
of the minority electrons. The form of the term linear in 
P2 is guarded by the rotational and time reversal symme- 
tries. What we need to prove is that the Taylor expansion 
indeed exists. This would be true if the dominant contri- 
bution to the vertex came from the "ultraviolet" region 
[in the same sense asused for derivation of Eq. (41)]. 

The only suspicious region where the non-analytic de- 
pendence of P2 can arise is the infrared integration in 
the diagrams which are reducible in one minority and one 
majority line, see Fig. 9. (a) and (b). Indeed, only in this 
case is the appearance of the overlaping poles possible, 
which may lead to the nonalyticity. Let us, however, ex- 
amine the expressions for both those diagrams in more 
detail. Let us write their analytic expression: 



Fig. 9. (a) + Fig. 9.(6) ~ j dilcfq x 

X V iq+,n+) V iq-,n_) G2 (e2 - O, p2 - q)) 
X [Gi (ei 17, pi q) Gi (ei - 17, Pi - q)] 

q±=q±|, n± = n±^ (44) 

where the screened potential is given by Eq. (15). The 
dangerous contribution may come only from the pole part 
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Pi+k/2 P2+k/2 
(a) 



I>,+ k/2 P2- k/2 

(b) 



FIG. 9: Lowest order diagrams for ri2. 



of the majority Green'S function (26). For \p\ = ppi we 
have 



Gi (e, p + q) = -Gi (-e, pi - q) : 



(45) 



which is nothing but the electron-hole symmetry for lin- 
earized spectrum. Therefore, the two terms in brackets 
in Eq. (44) cancel each other at ei, ^1 = 0, and therefore 
no infrared contribution is possible. All the corrections 



associated with the electron- hole asymmetry of major- 
ity electrons have at least one extra power of the Fermi 
energy epi in the denominator. 

It is clear that the same argument about the electron- 
hole symmetry remains valid even if the interaction lines 
V are replaced by the dressed vertices (43), and there- 
fore, the form (43) persists in all of the orders of the 
perturbation theory. 

So far, we established that the infrared part of the di- 
agrams involving mutual scattering of one minority and 
one majority electrons does not contribute because of 
the electron-hole symmetry of the majority system. It is 
clear that the above reasoning can be applied to higher 
order vertices as well. For further analysis, we will need 
only the 6-point vertex involving one majority and two 
minority electrons. Repeating all of the above consider- 
ation and taking into account the antysimmetricity (34) 
with respect to the permutations of the minority elec- 
trons, we find with logarithmic accuracy (the majority 
electrons are assumed to be on shell): 



^ir}^^ (P' 



out 
1 



■pout -pout, -pin pin pm\ 



b 



Fl 



(3) 1 

7i22 In 



PFl 



max(p2 



^in. ^out ^out^ 



V 



where 7^22 is the coefficient of the order of unity; we 
will not need its value in the subsequent consideration. 
This formula may be understood as the dependence of 
the prefactor in the two-particle vertex function (42) on 
the density of the majority electrons. 



5. Vertex functions W 



from Fig. 7. explicitly. We will see later that the value 
of W on shell is directly related to the Fermi liquid con- 
stants. We solve the diagramatic equation in Fig. 7. at 
Q "C Pfi, and iv <C Vpippi. Under this condition, we may 
replace Bj(P,Q) — > Sj(P,0) and use the Ward identity 
(37). Moreover, on majority shell i?i(P,0) = l/Zi and 
for e < VFiPFi and p < pfi, -B2(P,0) = I/Z2. This 
yields at u <^ vpiq 



The established dependence of the vertex functions 
r'lii ^22 enables us to find the dependence of functions W 



W||(Pi,P2;a;,q) 



1 + 



(0) 

11 



qvFi 



2^171^ 
Z2iyi 



2^1712 1 



^11 



+ (^1712) 



m 



The leading correction to Eq. (47) comes from Eq. (42) 
and it has the estimate (pi/p^i)^ ln(pf /p|^i). The mo- 
mentum dependence of the vertices B give a correction 
of the order of {pi/pri)^- We neglect the corrections of 
this kind. 

The other function we need in further calculations is 
^iK^i) P2; q). The Coulomb interaction line does 



not flip the spin of the electron and therefore this object 
is determined solely by the irreducible vertex (43). We 
find 



Wl2/p p.,, ^_ (Pl+P2)'-q' (2) f s 

1^12lPl)P2,W,q) - — -2 77i2 . (48) 

^1 1^1 \PFl\ 
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The minus sign here is associated with the change of the 
direction at which the irreducible vertex enters the dia- 
gram. 

In the following subsection we will use Eq. (47) to find 
the value of the Fermi liquid parameters. 



C. System near full polarization 




In this subsection we will take the following route. 
First, we calculate the changes in the Green function 

due to the finite density minority electrons but consid- 
ering their spectrum unchanged at n2 = 0. Second, we 
compute change in the vertex functions due to the finite 
density. Finally, we use the vertex functions to recalcu- 
late the spectrum of the minority and majority electrons, 
thus determining the Fermi liquid function. The small 
parameter justifying the procedure is 722/-^ -C 1. 



1. Green functions 

For nonzero minority electron density n2 > 0, n2 ^ 
rii we begin our discussion by considering the minority 
Green function G2 (hereinafter Green functions and n- 
point functions with tilde are understood at ^2 > 0, while 
the absence of the tilde implies n2 = 0). Finite density of 
majority electrons leads to the appearance of the positive 
chemical potential jj,2 and the shift of the pole in Eq. (29) 
to the upper semiplane, at Re e < 0. This change is 
described as 



62 (e, p) = G2 (e + /t2 , P) + ^G2 (e, p) + 60^ (e, p) • (49a) 

The second term in Eq. (49a) originates from the quasi- 
particle pole of the Green function 



6G2 = i27rZ2S (^e + il2- {pf2 - p) , (49b) 

where 6{x) is the Heaviside step-function. Here we ne- 
glected the correction to the parabolic spectrum and will 
restore this dependence later on. Within this approxima- 
tion, fl2 =p|,2/2"^2- 

The last term in Eq. (49a) gives zero while integrated 
over e only in the vicinity of the minority Fermi surface. 
However, this term is not an analytic function of e in the 
upper semiplane and therefore it gives finite contribution 
to the density. In the leading order in 722 it can be found 
from Fig. 10 and equals to 



SGr{e,p) = 

i [G2{e, p)]' J Til (P, Pi; P. Pi) '^^2(Pi)^ 



r 



FIG. 10: Contribution to 5E2. The internal line denotes 5G2- 



In this expression we implied that only p^ pf2 will con- 
tribute to the observable quantities and that is why we 
put two argument of the vertex function to zero. For the 
same reason, the chemical potential /i2 can be neglected 
in the argument of the Green function. 

The presence of the smooth term (49c) is crucial for the 
gauge invariance of the theory. In particular, it provides 
the cancellation of gauge-noninvariant factor Z2 from the 
observable quantities. As an example, we calculate the 
electron density n2. Using Eqs. (49) and the fact that 
G2{e) is analytic for Ime > 0, we find 



n2 = 



+m2^2 j (^r22 (P,0;P,0) G2(e,p)2 

_ ^2Pf2 



47r 



- n2 {Z2 - 1) , 



(50) 



where in the last transformation we used the Ward iden- 
tity (37), and 



de 



£^0 ^ Z2 ' 



It is seen from from Eq. (50) that the quasiparticle weight 
Z2 is cancelled and the Landau theorem 712 = p%2/^'^ 
holds. 



2. Vertex functions: Fermi liquid parameters and their 
corrections. 

To calculate the Fermi liquid constant Fjj (111112) one 
may start from the standard expression, see Fig. 11. 



Fij{n^2) = -Z,Z,v^W% (P; P; q; u; = 0) ; 



(51) 



n2^2r^^ (P,0;P,0) [G2(e,p)]^ 



(49c) 



P=l ^ — ,e = , q = pFjni -pirjn2. 



where the factor v\ is introduced in order to make the 

constants dimensionless. Notice, that with the vertex 
functions defined as in the previous section, the usual 
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R(niii2) = Vi Zi 



Pi 



Pi 



FIG. 11: Definition of Fermi- liquid constants. 



5G, 





on Fig. 12. (a)-(c). Contribution due to the irreducible 
6-point vertex is given by 



Fig. 12.(a) = 

j^SG^ (P) Till {P7\ Pr , P; PI" , P2", P) 



/ 



Using Eqs. (46) and (49b), we arrive to the estimate 



Fig. 12.(a) oc (na)^ , 



(a) 





and this contribution is neghgible. Contributions of 
Fig. 12. (b) and (c) are not small separately but their 
sum is. Presense of (5-function in Eq. (49b) guarantees 
the momentum and energy transfer to be small, therefore, 
the expansion (43) for the vertex functions is legitimate. 
We thus find 



FIG. 12: Leading contributions to F12. The filled hexagon 
denotes the irreducible vertex FiH. 

problem on non-commuting limits w ^ 0;q ^ does 
not arise at all. 

For the majority electrons Eq. (51) is just a formal def- 
inition which does not bring anything new. Actually, this 
definition was already used in the derivation of Eqs. (15) 
and (18). For the minority electrons, Eq. (51) actually 
allows one to calculate the Fermi liquid parameters. Sub- 
stituting Eq. (48) into Eq. (51) we find 



Fig. I2.(b)+Fig. 12.(c) = (^)7 ^^^G.iP) 
x[Gi(e,p + q) + Gi(-e,p-q)]. (53) 

Contribution from the pole parts of the majority Green 
functions is cancelled due to the electron-hole symmetry, 
compare with Eq. (44), and we obtain the finite contri- 
bution proportional to at least the first power of n2. 



4- Corrections to minoroty Fermi liquid parameter F22 



F12 (nUTa) = Z1Z2 7i2 + — ^12 "i • na . (52) 
V Pfi J 

In order to find F22 we substitute Eq. (47) into 

Eq. (51). Using zero angular harmonics of Eq. (52) to 

f 2) 

eliminate constant 7^2 , we obtain Eq. (22). 

The derivation presented here is still too cavalier. The 
expression (22), for instance, contains the dependence on 
the density of the majority electrons. However, this de- 
pendence came solely from the momentum dependence of 
the vertex functions W|| which was calculated at n2 = 0. 
Finite density ri2 leads to the contributions from the 
closed loops of the minority electrons and thus to modi- 
fication of To prove the legitimacy of retaining the 
second term in Eq. (22), one has to prove that the mod- 
ification of W|| is higher order in n^/N. We turn to the 
corresponding proof now. 

3. Corrections to minority-majority Fermi liquid parameter 

F12 

The correction to the Fermi-liquid constant F12 comes 
only from the irreducible vertex. This correction is shown 



Similarly, corrections to the Fermi liquid constant F22 
are determined by the six-point irreducible vertex, see 
Fig. 13. (a) and by the sum of six graphs determined by 
4-point vertices W2^2^ from Eq. (47), see Fig. 13. (b)-(g). 
Contribution of the first graph is estimated with the help 
of Eq. (40) for n = 3 and Eq. (49b): 

Fig. 13.(a) = 

/ ^^^^ (pr, pr*, p; Pi", Pi", p) 

/„in T-^^^^ fr-^out nOut\ 
^ \Pl ~ P2 j • IPl ~ P2 ) 

~ ^^ib^i]^ 

where the external momenta are assumed to be on shell. 
This result can be understood as the correction to the 
argument of the logarithm in Eq. (42). Because we al- 
ready neglected the term (42) in Eq. (47) as produc- 
ing (n2/A^) ln(r?,2/A^) correction, keeping the correction 
of Fig. 13. (a) would be also beyond the accuracy of the 
calculation and we can neglect it. 

Each term in the reducible graphs Fig. (b) - (e) is not 
small but their sum is. One finds 
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Fig. 13.(b) + • • • + Fig. 13.(c) = 
W22 [Pi + ^,P2 + §;Q 



(2n 



J [SG2{F2 + Q)G2(Pi + Q) + G2(P2 + Q)(5G2(Pi + Q)] 



W^22 Q + 



p. Pi 



-;r'2 



Pi 



1 2 



(54) 



The Green function SG2 from Eq. (49b) restrict the inte- 
gration in the infrared region, where one can use Eq. (47) 
for function W22, and Eqs. (29) and (30) for the Green 
function G2. The constant part of W22 immediately can- 
cels, and we repeat the dimensional analysis which lead 



us to estimate (40). This immediately yields the result 
proportional to n2, which can be neglected. 

The last two diagrams are analyzed in the same man- 
ner: 



Fig. 13.(f) +Fig. 13.(g) = J ^J(Pi;P2;Q) 



T(Pi; P2; Q) = [<5G2(P2 + Q)G2(Pi - Q) + G2(P2 + Q)SG2{Pi - Q)] W22 (Pi + ^, P2 + ^; Q 

/ Q Q \ /P1 + P2 P1 + P2 

M^22 fPi + ^,P2 + -^;Q] -W^22 ( 2 2 + 

Shifting the integration variable, we replace I(P 1 ; P2; Q) [I(Pi; P2; Q)-|-I(Pi; P2; Q-I-P2 — Pi)]/2 in the integrand 
and find 



Fig. 13.(f) + Fig. 13.(g) = i 1 ^ [JG2(P2 + Q)G2(Pi - Q) + G2(P2 + Q)<5G2(Pi - Q)] 

/ Q Q \ /P1 + P2 P1+P2 
W^22 fPi + ^,P2 + ^;Qj -W^22 ( I -Q, I ' +Q;Q + Pi-P2 



Similarly to Eq. (54), it yields the result proportional to 
the first power of n2. 



5. Final remarks 

Let us now discuss the effect of the finite density of 
minority electrons on the majority electrons. Firstly, one 
observes that the minority polarization operator is not 
small: 112 = 1^2 at momentum transfer q < 2pp2 and 
decays as 112 — v2{pf2/q)'^ larger momenta. Naive 
calculation of the correction to the majority Fermi liquid 
parameter Fii{9) gives large contribution to the forward 
scattering for angles < Pf2/pfi- This contribution, 
however, is cancelled out from any closed loop for the 
majority electrons, see Fig. 14. The easiest way to see it, 
is by noticing that any closed loop with arbitrary number 
of scalar vertices vanishes due to the gauge invariance if 
one of the momenta equals to zero. Therefore, one may 
gauge away those contributions from the very beginning, 
and leave only processes with momentum transfer of the 



order of ppi . In this region, correction to the majority 
loops due to the presence of a nonvanishing density of 

minority electrons is already small as [pf2/pfi]'^ , which 
ones again results in the correction proportional to the 
first power of n2- 

Before concluding this subsection wc have to derive 
Eq. (8b) governing the renormalization of the minority 
mass. Taking the expression of F12 into account from 
Eq. (52) wc sec that the large ratio Pfi/pf2 in Eq. (24) 
is cancelled in the angle-dependent term, and hence we 
obtain a renormalization of the inverse mass 1/m by a 
quantity of order 1/mi. Utilizing the expression (43) of 
ri2, we can immediately see, that further corrections to 
the first harmonic of are smaller by at least 712 /N and 
hence result in a nonsingular renormalization of the mass. 
On the other h and, th e first harmonic of F22 is nonana- 
lytic, of order \fn2jN . This implies that at zero minority 
density the renormalization of the minority mass is due 
to the majority "background" described by F12, and at 
finite densities it can be attributed to the residual inter- 
action (described by F22 via the real part of V) between 
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(e) (f) (g) 



FIG. 13: Leading contributions to F22. The filled hexagon 
denotes the irreducible vertex T222 ■ For each of the diagrams 
(b)-(g) there exists a counterpart obtained by swapping 5G2 
and G2. 




1 1 



FIG. 14: Example for the cancellation of the singularity in 
forward scattering. The majority loop can be further dressed 
in all possible ways. 



minority quasiparticles. 

To summarize the content of this subsection, we 
showed by expUcit analysis of all orders of perturbation 
theory, that the naive derivation of the Fermi liquid con- 
stants for the two fluid model, Eq. (22), and the minority 
DOS, Eq. (8b), is parametrically justified at n2 4C N. 
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APPENDIX A 

In this Appendix we justify the form of the minority 
self-energy given in Eq. (31) and calculate the dimen- 

sionlcss function f{x). Clearly, "ultraviolet" self-energy 
diagrams [in the sense discussed after Eq. (40)] can be 
expanded in Taylor series in both and e. These di- 
agrams are responsible for the renormalization of the 
chemical potential (to some negative or zero value), the 
minority mass (which is assumed to be positive) and the 
quasiparticle weight Z2- All nonanalytic behavior of the 
zero minority density self-energy must come from the "in- 
frared" , where a perturbation expansion becomes possi- 
ble in terms of the small nonconstant part of the screened 
interaction. 

Our strategy is the following. We first calculate the 
contribution to the imaginary part of the minority self- 
energy S2(e, p) coming from the diagram in Fig. 15. 
Then, we will argue that all other infrared diagrams are 
smaller by a factor of at least [p/pfi]^- The contribution 
of Fig. 15. to the self-energy is 

Me,p)=i J / ^<5y(a;,q)G2(e + a;,p + q)B|. 

For the imaginary part we obtain 



IV. SUMMARY 

In conclusion we studied the almost fully polarized 
two-dimensional electron liquid. We have shown that as- 
suming (i) the stability of the majority Fermi-liquid and 
(ii) the positive renormalization of the minority mass, 
no matter how small the density of the minority spin 
electrons is, the Fermi two-liquid description is always 
consistant. Moreover, microscopic analysis made it pos- 
sible to find the connection between different Fermi liquid 
paramaters, thereby reducing the number of independent 
parameters by one. 

The established Fermi liquid description enable us to 

predict important relations between different thermody- 
namic observables of the system, see Subsec. II B. 



The theta function restricts the domain of integration to 
momenta (p-|-q)^ < 2m2e. Since e <C epi we can further 
simplify the above formula by (i) using the Ward identity 
Eq. (37), thus effectively setting B2 = 1/^2, and (ii) by 
expanding the screened potential V according to Eq. (38) 
for small energies and momenta and fov w <^ vpiq- This 
way we arrive to 

-I„.E,(,|p|) =^-L.J^. (A2) 

VFi\q.\ 
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FIG. 15: Lowest order contribution to the minority self- 
energy. The dashed line denotes the screened Coulomb inter- 
action 1/, see text, the solid line is the minority propagator 
Gi. Black triangles will be explained in Subsec. Ill A. 



Here the energies of the external legs were put to zero for 
simplicity. 

In the infrared region (g < Q) we deform the fre- 
quency integration contour, expand the renormahzed in- 
teraction y for w < vf\(1 ^ Vf\ according to Eq. (38), 
and introduce the dimensionless integration variables 
X = q/Q, y = — 2m2w/Q^. We obtain 



Calculating the integral and retrieving the real part of 
the self-energy in a way to make it analytic on the upper 
half-plane, we arrive to Eq. (31) 



-E2(e, IpI) = / f 



Pfi ' \2m2e 
with the dimensionless function / being 



fix) = ^ {5(1 - x)K{x'^^) - 2(1 + Ax)E{x'/^)} , 

(A3) 

where K{y) [E{y)] denotes the usual complete elliptic 
integral of the first [second] kind. 

Note, however, that the above calculated expression 
of the self-energy does not contain the full real part, 
sinc:c; an arbitrary analytic function of energy (with pos- 
sible momentum dependence) could still be added. This 
means that one cannot retrieve the real part of the 
sclf-encrgy, an inherently ultravialet quantity, from the 
infrared-related imaginary part. 

The contribution of higher order infrared diagrams is 
of order ma.x{e/epi, [p/ppi]'^} because of the smallness 
of the nonconstant part of the screened interaction V in 
the infrared region of momenta and energies. 



r2,2 

^ 2,2 



, Pi P2 



1 



V2 [PFl 



(B2) 



/ xdx I 
Jo Jo 



^PFl/Q 



dy 



y 



(y + x^)^ 



where the smooth, dimensionless real function 71 de- 
scribes the dependence on the external momenta. This 
integral is logarithmicaly divergent at the upper limit of 
the integration over y. Corrections may come from the 
further expansion of the potential V: these, however, 
result in a convergent integral, and arc therefore small 
compared to the leading logarithmic behavior. 

In the ultraviolet region {q > Q) the integration mo- 
mentum is always bigger than the external momenta, so 
we can Taylor expand the integral around pi = P2 = 0. 
Since the square bracket in Eq. (Bl) vanishes for zero 
external momenta, and because of rotational invariance, 

the expansion must start with Q^72 , 72 being 

another real dimensionless function. Now we can deform 
the contour of the frequency integration and expand V 
again. Introducing the dimensionless parameters x and 
y we obtain 



APPENDIX B 

In Appendix B we outline the derivation of Eq. (42), 
which was obtained using a simple dimension-counting 
argument. As already mentioned, the leading order con- 
tribution to the vertex r2'2 comes from the irreducible 
tree level diagrams Fig. 4. (a), (b), (d) and (e). 

Let us first pick diagram (a) and its crossed coun- 
terpart, (d) and write their contribution to the vertex 
r2|2(Pi,P2;Pi,P2) up to an overall factor as 

j (fq j duj G2(a;,q)V'(a;,q-pi) x (Bl) 

[G2(a;,q) y(a;,q-p2)- 

G2{uJ, q + P2 - T>i)V{u}, q - pi)] . 



T.2,2 . /Pi P2\ 1 

['^ f~^^^ { 1 \ 



an integral also logaritmically divergent at the upper 
limit. 

The other two diagrams Fig. 4. (b) and (e) give similar 
contributions. The dressing of the scalar vertices ensures 
the cancellation of the quasiparticle weight Zi of the mi- 
nority Green function (omitted in the above estimates for 
clarity). Consequently, one obtains Eq. (42). 
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